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acceleration. 
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1. Introduction 

The simplest homogeneous and isotropic cosmological model within the framework of 
general relativity is parameterized by two numbers: the Hubble constant Hq and the 
density parameter Qm, interpreted respectively as the average expansion rate and 
the average density of non-relativistic matter in the present universe. This model was 
generally considered a good description of the universe on the largest scales until the 
high redshift supernova observations in the late 90's [1]. With the latest data from 
supernovae [2,3], galaxy distributions [4] and anisotropics of the cosmic microwave 
background [5] the simplest homogenous and isotropic model would now lead to a 
highly contradictory picture of the universe, as the best fit values of this model for 
the average matter density demonstrate: 

• Cosmic microwave background: Qm ~ 1 
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• Galaxy surveys: VIm ~ 0.3 

• Type la supernovae: fi^ ~ 

The natural conclusion is that at least one of the assumptions of the model must be 
false. 

As is well known, the problem has conventionally been remedied by introducing 
the cosmological constant A or vacuum energy giving rise to an accelerated 
expansion of the universe. Indeed, the analysis of the cosmological data with both 
the vacuum energy and matter components included yields a consistent picture of 
the universe, known as the concordance ACDM model, with the following best fit 
values for the density parameters: 

• Cosmic microwave background: VLm + ~ 1 

• Galaxy surveys: Qm ~ 0.3 

• Type la supernovae 1 : Vl M ~ 0.3 and ~ 0.7 

Moreover, the data seems to require Q\ > at a high confidence level [2, 3]. 

Although the cosmological concordance ACDM model fits all the observations 
well, it is plagued by theoretical problems [6]. We do not have a theory that would 
explain, not to mention predict, the required value of the cosmological constant. 
Various dark energy models, which have been studied intensively in order to provide 
a dynamical explanation for the cosmological constant, are not compelling from the 
particle physics point of view and often require fine-tuning [6, 7]. Modifications of the 
general theory of relativity on cosmological scales appear to suffer from analogous 
problems; in fact, it has recently been argued that the cosmological constant seems 
to be essentially the only modification that fits all the cosmological data [8]. 

As a rejection of either matter domination or Einstein gravity leads to trouble, 
it is well motivated to study the validity of the third main assumption, the per- 
fect homogeneity. Undoubtedly, at least on small scales, large inhomogeneities exist. 
Although the potential cosmological consequences of the inhomogeneities were rec- 
ognized already around the same time when the homogeneous and isotropic models 
of the universe were first studied, their impact on the global dynamics of the uni- 
verse is still unknown (see e.g. [9]). Indeed, the role of the inhomogeneities has often 
been debated whenever there have been some ambiguities with cosmological obser- 
vations, for example, in the 90's when in the determination of the age of the universe 
there was a discrepancy between the implications of cosmological and astronomical 
observations [10,11]. 

1 In the supernova data analysis spatial flatness (<S=4> SIm + Oa = 1) has been assumed. However, 
these values are within the la error of the best fit values without the flatness constraint [2]: Q M <~ 
0.5 and ~ f . 
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Most recently, inhomogeneities have been invoked as the culprit for the apparent 
acceleration of the expansion of the universe 2 , in particular by virtue of their so- 
called backreaction on the metric (for a recent discussion on the issues involved and 
a comprehensive list of references, see [13]). The problems in these approaches have 
been mainly of technical nature, as reducing the symmetries of the metric rapidly 
complicates the calculations. A lesser difficulty arises from the fact that different 
inhomogeneities can lead to identical observations [14], so that not even ideal ob- 
servations of light could determine the inhomogeneities uniquely. As a consequence, 
there have been two conceptually distinct ways to approach the inhomogeneities: the 
first one examines the effect of the inhomogeneities on the expansion of the universe, 
whereas the second tries to determine their impact on the observations directly. 

In the first approach an effective description of the universe is constructed by 
averaging out the spatial degrees of freedom, i.e. the inhomogeneities [13,15-17]. 
As a result, one obtains averaged, effective Einstein equations which, in addition to 
terms found in the usual homogeneous case, include new terms that represent the 
effect of the inhomogeneities, called backreaction in this context. It has been argued 
that the backreaction might account for the accelerated expansion [13, 16, 18]. 

However, since we can only observe the redshift and energy flux of light arriving 
from a given source, not the expansion rate and the matter density of the universe 
nor their averages, one may wonder how the actual observables are related to the 
averaged equations. To wit, since we do not observe the expansion of the universe 
directly, its acceleration is also an indirect conclusion, arising from the fact that in 
the perfectly homogeneous cosmological models dark energy is required for a good fit. 
Consequently, there is no a priori reason to assume that the acceleration would be 
needed in the more general inhomogeneous models of the universe. Moreover, it is 
well recognized that in general it is not correct to integrate out constrained degrees 
of freedom as if they were independent. Indeed, the fact that we can make cosmo- 
logical observations only along our past light cone makes the observable universe a 
constrained system. 

The second approach avoids these problems by studying the effect of the inhomo- 
geneities directly on the observable light [14, 19]. This would be virtually impossible 
in the presence of generic inhomogeneities but can be done in some simpler mod- 
els, such as the spherically symmetric Lemaitre-Tolman-Bondi (LTB) model [20-22] . 
This model works best in describing smooth inhomogeneities at scales of 100 Mpc 
and larger; the spherical symmetry prevents to use it as a model for the random 
small scale lumpiness caused by galaxies, as noted in [23]. Indeed, the LTB model 
has been used to study the effect of the smooth inhomogeneities on the cosmological 
observations by several authors [23-32]; a common conception is that these models 
inevitably contradict the observed homogeneity of the large scale galaxy distribution. 

2 Inhomogeneities as an alternative to dark energy were first discussed in [12]. 
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Although spherical symmetry is probably an unrealistic assumption for the entire 
universe, the LTB model can be regarded as describing observations that have been 
averaged 3 over the celestial sphere and it is therefore useful at least on two counts. 
First, it serves as a simple testing ground for the effect of the inhomogeneities on the 
cosmological observations. Second, as the fits can be performed unambiguously, it 
can be used to study the connection of the backreaction driven effective acceleration 
to the observations; both of these points will be examined in this work. 

As mentioned in the beginning, two independent parameters (Hq and Qm) are 
needed to uniquely define the homogeneous matter dominated universe. In the pres- 
ence of inhomogeneities, the values of these two quantities are needed at every spatial 
point; that is, the inhomogeneous dust models are defined by two functions of the 
spatial coordinates: H (x l ) and VLm{x 1 ). As a consequence, there are inhomogeneities 
of two physically different kind: inhomogeneities in the matter distribution, and in- 
homogeneities in the expansion rate. Although their dynamics are coupled via the 
Einstein equation, as boundary conditions they are independent. This opens up the 
possibility for a universe with inhomogeneous expansion but homogeneous present- 
day matter distribution; a model of this kind could potentially fit the supernova data 
as well as the galaxy surveys without invoking dark energy. 

The paper is organized as follows. Sect. @ discusses the general properties of the 
spherically symmetric, inhomogeneous LTB model. In Sect. [|we fit the LTB model 
with inhomogeneous expansion but homogeneous present-day matter distribution 
to the supernova observations. For completeness, we consider both the cases of 
pure dust and dust plus the cosmological constant. The initial conditions of these 
models are discussed in Sects. [T7] and pT8| . In Sect. [| we evaluate the expansion 
rate and shear in the appropriate averaged (the so-called Buchert [16]) equations and 
demonstrate that an accelerated average expansion is not needed to fit the supernova 
data. Finally, Sect. [| contains our conclusions. 

2. Spherically symmetric inhomogeneous LTB model 

Since Einstein's equations form a set of six independent non-linear second-order 
partial differential equations, it is impossible to treat the universe exactly with com- 
pletely generic inhomogeneities. On the other hand, the perfectly homogeneous FRW 
model does not fit the observations without a fine-tuned cosmological constant or 
some other modification. Therefore, we make the next simplest approximation and 
use the spherically symmetric but inhomogeneous LTB model. The aim is to extract 
the leading order effects of the inhomogeneities on the supernova observations as 
well as to demonstrate the pitfalls one may encounter when using averaged Einstein 
equations. 

3 This is naturally only a crude approximation; a more correct way would be to first solve the 
observables in a nonspherical model and only then average out their angular dependence. 
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The LTB model has been used for the supernova data fitting several times before 
[23,24,26-30]. However, there is a crucial physical difference between these models 
and ours. Namely, the earlier works have had inhomogeneities in the present-day 
matter distribution whereas we will focus on models with inhomogeneities only in the 
expansion rate. Moreover, to ease the comparison between the FRW model, familiar 
to all cosmologists, and the less-known LTB model, we will rewrite the equations in a 
form where the connection of the physically equivalent quantities between these two 
models becomes very transparent. To introduce the new notation, we will rederive 



the main results of the LTB model in Sects. |2.1| and |272 . 
2.1 The Lemaitre-Tolman-Bondi metric 

Let us consider a spherically symmetric dust universe with radial inhomogeneities 
as seen from our location at the center. Choosing spatial coordinates to comove 
{dx l /dt = 0) with the matter, the spatial origin (x* = 0) as the symmetry center, 
and the time coordinate (x° = t) to measure the proper time of the comoving fluid, 
the line element takes the form 

ds 2 = -dt 2 + ^ dr 2 + A 2 (r, t) (d6 2 + sin 2 Bdy 2 ) , (2.1) 
1 — k(r) 

where k(r) is a function associated with the curvature of t = const, hypersurfaces, 
the scale function A(r,t) has both temporal and spatial dependence, and we use 
the following shorthand notations for the partial derivatives: ' = and = J^. 
This metric was first studied by Lemaitre [20], Tolman [21] and Bondi [22]; later, it 
has been used in various astronomical and cosmological contexts [9]. Note that the 
homogeneous and isotropic FRW metric is a special case of Eq. (|2.1|) , obtained in 
the limit: A(r,t) —>■ a(t)r and k(r) — > kr 2 , where a(t) is the FRW scale factor and k 
is the curvature constant. 

The energy momentum tensor in the above defined coordinates is given by 

T^ = -p M {r,t)6Z5° u -p A 5» u , (2.2) 

where Pm{ t i t) is the matter density, w M = Sq represents the components of the 4- 
velocity-field of the fluid and we have kept the vacuum energy p A for generality. 
Although the fluid is staying at fixed spatial coordinates, it can move physically in 
the radial direction; this movement is encoded in ^/gn = A'(r, t) / yT — k(r). 

When Eqs. ( J2.1| ) and ( |2.2j ) are applied to the Einstein equation, G li v = SttGT^, 
two independent differential equations arise: 

A 2 + k(r) 2AA' + k'(r) , . 

A t + J^ 1 = ^G(p M + p A ) (2.3) 

A 2 + 2AA + k{r) = 8nGp A A 2 . (2.4) 
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The first integral of Eq. ( [2.4|) is 



A 2 



F(r) 8nG k(r) 
V ^^Pa 



(2.5) 



A 3 ' 3 riv A 2 

where F(r) is a non-negative function. Substituting Eq. (|2.5|) into Eq. ( |2.3| ) gives 

F' 



A 1 A 2 



8nGp 



M 



2.6) 



By combining Eqs. 
equation 



and (|2.4| ) we can construct the generalized acceleration 



2 A I A' 

3 A + 3 A 7 



4ttG 



(2.7) 



This equation tells that the total acceleration, represented by the left hand side, is 
negative everywhere unless the vacuum energy is large enough: p A > pm/2- However, 
it does not exclude the possibility of having radial acceleration (A'(r, t) > 0), even in 
the pure dust universe, if the angular scale factor A(r, t) is decelerating enough and 
vice versa. Already a simple example like this demonstrates how the very notion of 
the acceleration becomes ambiguous in the presence of the inhomogeneities [33]. 

The boundary condition functions F(r) and k(r) are specified by the exact nature 
of the inhomogeneities. Their relation to the more familiar quantities — the Hubble 
constant Ho and the density parameter Om — can be recognized by comparing Eq. 
(|2.5|) with the Einstein equation of the homogeneous FRW model 



H 2 (t) = 



a 2 {t) 
^(t) 



Hi 



n 



M 



OoV 

a J 



(1-fi/ 



n 



M 



OoV 

a J 



X8) 



where a = a(t ). Thus, the comparison of Eqs. 
define the local Hubble rate 

A(r,t) 



H(r,t) = 



A(r,t) 



31) and (E.8I) motivates us to 



(2.9) 



and local matter density through 



F(r) 



H 2 (r)n M (r)Al(r) 



(2.10) 



with 



k(r) = H*(r)(Q M (r) + fl A (r) - l)A 2 (r) , (2.11) 

where A (r) = A(r,t ), H (r) = H(r,t ), and n A (r) = 8nGp A /3H 2 (r). With these 
definitions Eq. (|2.5| ) takes the physically more transparent form 



H\r,t) = Hl{r) 



(2.12) 
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where O c (r) = 1—Q\(r)—Q^j(r). The difference between the conventional Friedmann 
equation (|2.8| ) and its LTB generalization, Eq. (|2.12|) , is that all the quantities in 



the LTB case depend on the r-coordinate. This is true even for the gauge freedom 
of the scale function: In the FRW case the present value of the scale factor a(to) can 
be chosen to be any positive number. Similarly, the corresponding present-day scale 
function A(r,t ) of the LTB model can be chosen to be any smooth and invertible 
positive function. For the rest of this paper, we will choose the conventional gauge 

A(r,t )=r. (2.13) 

Although the vacuum energy density p\ is constant, its value in the units of 
critical density fl\(r) = p A / p cr it( r ) is not. This is because the critical density itself 
has spatial dependence: p CT it(^) = 3i/o(r)/87rG. The converse is also true: if e.g. 
^m(^) = constant, the matter distribution pu itself has spatial dependence as long 
as H Q {r) 7^ constant. 

Integrating Eq. (|2.12| ) gives 

If 1 dx , 

t -t= — — / = . 2.14 

H (r) Jm0l x^n M (r)x- 3 + fi c (r)x" 2 + fi A (r) 

For any space-time point with coordinates (t,r,9,<p), Eq. fl2.14| ) determines the 
function A(r, t) and all its derivatives. Thus the metric in Eq. Q2.ip is specified and 
given the inhomogeneities, all the observable quantities can be computed. 

2.2 Relation of the inhomogeneities to the observations of light 

To compare the inhomogeneous LTB model with the supernova observations, we 
need an equation that relates the redshift and energy flux of light with the exact 
nature of the inhomogeneities. For this, we must study light propagation in the LTB 
universe. We will again rederive the appropriate equations for notational clarity; a 
more general derivation for an off-center observer can be found in [34]. 

From the symmetry of the situation, it is clear that light can travel radially, that 
is, there exist geodesies with d6 = dip = 0. Moreover, since light always travels along 
null geodesies, we have ds 2 = 0. Inserting these conditions into the equation for the 
line element fl2.1|) , we obtain the constraint equation for light rays 

g f ?™ (2.i5) 

du du yT - k(r) 

where u is a curve parameter and the minus sign indicates that we are studying 
radially incoming light rays. 

Consider two light rays with solutions to Eq. ( |2.15| ) given by t\ = t(u) and 
h = t(u) + X(u). Inserting these to Eq. ( |2.15| ) we obtain 

d dt(u) dr A'(r,t) 

-Hi = — = j , 2.16 

du du du y/i — k(r) 
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= ^M. + rfA (^) = + dA M ( 2 17) 

du du du du a/1 — k(r) du 

drA'(r,*(u) + A(u)) rfr A'(r, t) + A'(r, t)X(u) 

1~^2 = ; = = : = , (2.18) 

du du y/i - k(r) du ^1 - k(r) 

where Taylor expansion has been used in the last step and only terms linear in A (it) 
have been kept. Combining the right hand sides of Eqs. fl2.17| ) and ( |2.18| ) gives the 
equality 

d\(u) = dr A'(r,t)X(u) ^ 
du du y/l - k(r) 

Differentiating the definition of the redshift, z = (A(0) — X(u))/X(u), we obtain 

dz _ dX(u) A(0) _ dr (l + z)A'(r,t) 
du du X 2 (u) du yj\ — k(r) 

where in the last step we have used Eq. ( 2.19| ) and the definition of the redshift. 
Finally, we can combine Eqs. (|2.15|) and (|2.20|) to obtain the pair of differential 

equations 

dt -A'{r,t) 



dz (l + z)A'{r,t) 
dr _ y/1 + Hl(r){\ - Q M (r) - Q A (r))A 2 (r) 



(2.21) 
(2 22) 

dz (l + z)A'(r,t) ' V ' ' 

determining the relations between the coordinates and the observable redshift: t(z) 
and r(z). 

Now that we have related the redshift to the inhomogeneities, we still need 
the relation between the redshift and the energy flux F, or the luminosity- distance, 
defined as d^ = ^L/AnF, where L is the total power radiated by the source. The 
desired relation is given by [35] 

d L (z) = (l + z) 2 A(r(z),t(z)) . (2.23) 

As the relations t(z) and r(z) are determined by Eqs. Q2.21 ) and ( |2.22| ) and the scale 
function A(r, t) by Eq. ( |2.14|) , using Eq. ( p.23|) one can calculate di for a given z. 
All of these relations have a manifest dependence on the inhomogeneities (i.e. on the 
functions H (r) and flM{ r ))- What remains is a comparison of Eq. ( |2.23| ) with the 
observed d^z). 

In the FRW model the parameters that best describe our universe are found by 
maximizing the likelihood function e -x 2 {H ,Q, M ,o, A ) cons tructed from the observations. 
However, to find the boundary conditions of the LTB universe that best describe our 
universe, we should in principle maximize the likelihood functional e~ x2 ^ Ho ^' nAl( - r ^ . 
In practice, this is impossible. Therefore we will consider some physically motivated 
types for the functions Ho(r) and VIm{t) that contain free parameters; these are then 
fitted to the supernova observations by maximizing the leftover likelihood function. 
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3. Models with inhomogeneous expansion 

In this Section, we will discuss four models where the form of the boundary condition 
functions H (r) and VLm{t) has been specified and fit them to the data of the Riess 
et. al. gold sample of 157 supernovae 4 [2]. An extensive cosmological data analysis 
would also take into account the galaxy surveys and the anisotropies of the cosmic 
microwave background. However, we will be satisfied to give only qualitative argu- 
ments as to why our models could have the potential to fit these data sets as well; 
discussion of the potential fit for the CMB data will be given in Sects. and §. 

According to the observations, galaxies seem to be rather evenly distributed in 
the universe [4]. Indeed, the usual objection to the inhomogeneous models is that 
they contradict the observed homogeneity of the large scale structure [28]. However, 
as stressed in Sect. 0, the matter distribution of the present universe can be homoge- 
neous even though the expansion rate would have spatial variations. This motivates 
us to consider models with a uniform present-day matter distribution. 

When modelling the universe as perfectly homogeneous, its expansion has to 
accelerate in order to fit the supernova data. Mathematically, acceleration means 
that the second time derivative of the scale function is positive, but as can be seen 
from Eq. ( |2.7| ), this is not possible without the cosmological constant or some other 
form of dark energy. But then again, the observations are made along the past light 
cone, and what affects the observations is the variation of the dynamical quantities 
along the past light cone, not just the time variation. This is naturally true in the 
homogeneous universe as well, but as the time variation differs from the variation 
along the light cone only in the presence of inhomogeneities, one does not usually 
bother to make the distinction. However, here the difference is essential. 

The directional derivative along the past light cone reads as 



where the approximation in the last step is more accurate for the small values of r, 
but is qualitatively correct even for larger r. The main content of Eq. ( |3.1|) is that 
from the observational point of view, the negative r-derivative roughly corresponds 
to the positive time derivative. This is natural since by looking at a source, we 
simultaneously look into the past (i.e. along the negative t-axis) and spatially further 
(i.e. along the positive r-axis). So to mimic the acceleration, i.e. for the expansion 
rate to look like it would increase towards us along the past light cone, the expansion 
H Q (r) must decrease as r grows: H' (r) < 0. 

4 Note that we have chosen the inner luminosity, or the total radiation power, of type la super- 
novae as in [2]. 



d d dr d d \Jl — k{r) Odd 
dt dt dt dr dt A'(r,t) dr dt dr ' 
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With the above given arguments in mind, we have chosen the following form for 
the boundary condition functions: 



H (r) = H + AHe~ r/r ° , 

^m(t") = Qq = constant , (3.2) 

where H, AH, r and Q are free parameters determined by the supernova obser- 
vations and the exponential has been chosen for simplicity. In the model of Sect, 
all the four parameters are left free whereas in Sect. p.2| we have fixed Oq = 1. 



For generality, the cosmological constant has been included in Sect. |3l| for compu- 
tational simplicity and to facilitate the comparison with the model of Sect. |37| we 
have also set: Ojw(r) + Q\(r) = 1. 

As explained in Sect. [2.1| , the present-day matter distribution in the models 
with Qm{ t ) — const, is not perfectly uniform since the critical density depends on 
Hq{t). Hence we have also studied a model with pu^^o) = const, in Sect. |3T1[ 
The discussion of all the fits is given in Sect. [3.6| and the possible problems with the 



initial conditions of these models are then discussed in Sects. |3T7| and 



3.1 Inhomogeneous expansion and dust with O^W = const. 

For fixed values of the parameters (H, AH, r , Q ), one can compute the luminosity- 
distance-redshift relation of Eq. ( j2.23| ). By repeating this computation for different 
values of the parameters we search for the maximum value of the likelihood function 
e -x 2 (HAH,r ,n ) j_ e _ ^he minimum value of 



x 2 = E { dfB{Zn) " dlj{Zn> ) > ( 3 - 3 ) 



n=l 



where d°^ s {z n ) is the observed luminosity-distance for a source with redshift z n and 
a n is the estimated error of the measured d°^ s . In this way, we find that the best fit 
values for the parameters in this model are: 

• H + AH = 66.8 km/s/Mpc 

• AH = 10.5 km/s/Mpc 

• r = 500 Mpc 

• tt = 0.45 

• Goodness of the fit: x 2 = 172.6, x 2 /157 = 1.10 

The confidence level contours with Q and H fixed to their best fit values are shown 
in Fig. |l|. For comparison with the homogeneous case, the best fit nonflat ACDM 
has tt M = 0.5, tt A = 1.0 and X 2 = 175, x 2 /157 = 1.11. 
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2600 



Length scale of the variation in Mpc (parameter r Q ) 



Figure 1: Confidence level contours in the LTB model with 0,M( r ) = constant = 0.45 and 
H (r) = 56.3 km/s/Mpc + AHe~ r / r °. 



3.2 Inhomogeneous expansion and dust with ^A/(r) = 1 

The boundary condition functions have here the same form as in the model of Sect. 
|3.1| , but this time we have fixed VLm{t) = 1 in analogy with the flat FRW model. A 
great virtue of this model is that the scale function can be explicitly solved from Eq. 
Q2.14| ). In this case, the best fit values for the parameters are: 



• H + AH = 65.5 km/s/Mpc 

• AH = 16.8 km/s/Mpc 

• r = 1400 Mpc 

• Goodness of the fit: x 2 = 176.3, % 2 /157 = 1.12 

The confidence level contours with H fixed to its best fit value are displayed in Fig. 
|2|. For comparison with the homogeneous case, the best fit flat concordance ACDM 
model has Q M = 0.3, fi A = 0.7 and x 2 = 177, x 2 /157 = 1.13. 
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Figure 2: Confidence level contours in the LTB model with Q^(r) = constant = 1 and 
H (r) = 48.7 km/s/Mpc + AHe~ r / r °. 

3.3 Inhomogeneous expansion with ^m(t) + ^A(r) = 1 

We have now seen that inhomogeneous models of the universe can fit the supernova 
data without the cosmological constant. On the other hand, the homogeneous model 
with the cosmological constant fits the data as well. Thus, it is natural to ask whether 
an even better fit could be obtained with both the cosmological constant and the 
inhomogeneities. Let us therefore require that ^m(^) + ^a(V) = 1 with the same 
form for the function H (r) as before. Then we find 

• H + AH = 66 km/s/Mpc 

• AH = 8 km/s/Mpc 

• r = 600 Mpc 

• = PA/Pcrit = 0.33 for fixed p crit = 3H 2 / (8nG) with H = 66 km/s/Mpc 

• Goodness of the fit: X 2 = 173.4, x 2 /157 = 1.10 

Comparing this model with the homogeneous flat ACDM (x 2 /157 = 1.13) and the 
analogous inhomogeneous dust model of Sect. [T2] (x 2 /157 = 1.12) suggests that the 
mixture of the inhomogeneous expansion and the vacuum energy does not give a 
better fit than either of them separately. 
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Figure 3: Confidence level contours in the LTB model with perfectly uniform present-day 
matter density: H (r) = 57 km/s/Mpc + AHe~ r / r \ tt M (r) = 0.29(67 km/s/Mpc) 2 /#o (r). 

3.4 Inhomogeneous expansion and dust with pu{ r i to) — const. 

Let us then consider a model with strictly uniform present-day matter distribution: 
pM^ito) = constant. As can be seen in Eq. ( |3.8| ), the condition for the uniformity 
is: Hi {r)VL M {r) = constant. Therefore, we choose the boundary condition functions 
to be of the form 



The data analysis then gives the following best fit values: 

• H + AH = 67 km/s/Mpc 

• AH = 10 km/s/Mpc 

• r = 450 Mpc 

• n = 0.29 

• Goodness of the fit: X 2 = 172.6, x 2 /157 = 1.10 

The confidence level contours with Q and H fixed to their best fit values are displayed 



H (r) = H + AHe~ r/ro 




(3.4) 



in Fig. 0. 
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Figure 4: The functions Om(t*) (red line), Hq(t) in the units of 100 km/s/Mpc (blue line) 
and the physical matter density today Pm(^^o) i n the units of the critical density at the 
origin p c rit(0) (green line) for the best fit values in the model with simultaneous Big Bang. 



3.5 Inhomogeneous expansion with simultaneous Big Bang 



Finally, let us consider a model with simultaneous Big Bang, i.e. a spatially constant 
age of the universe ( |A.1| ). This constraint leaves us with only one free function. 
Hence, to maintain the radially decreasing Hubble expansion, we choose the bound- 
ary condition functions to have the form 



H (r) = H 



'1 - n M (r) - fi A/ ( r )arsinh A /^g^ 



(l-0 M (r))^ 



(l + 5e- r / r °) 2 



(3.5) 



The data analysis gives that the best fit values for these parameters are: H = 76.5 
km/s/Mpc, 5 = 1.21, r = 1000 Mpc, Q = 0.92, goodness of the fit: \ 2 = 175.5, 
X 2 /157 = 1.12. These values imply that the Hubble function Hq(t) varies from 
the value Ho(0) = 65 km/s/Mpc near us to its asymptotic value H (r ^> r ) = 52 
km/s/Mpc, as shown in Fig H. Eq. (|A . 1|) then gives the age of the universe simply 



as t age = 1/H = 12.8 Gyr. The values are similar to the model of Ref. [26]. 
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3.6 Discussion of the fits 



As proved by Mustapha, Hellaby and Ellis, any isotropic set of observations can 
be explained by appropriate inhomogeneities in the LTB model [14]. Hence, the 
fact that we could find boundary conditions fitting the supernova observations is not 
surprising. Instead, the novel aspect here is the form of these functions that gives the 
best fit, as well as their physical interpretation: inhomogeneities in the expansion rate 
but homogeneous present-day matter density. The particular interest lies in the fact 
that the analysis is not only in qualitative agreement with the observed homogeneity 
in galaxy surveys, but also gives a similar value for the present-day matter density: 
Qm ^ 0.4. Of course, the matter power spectrum should be reanalyzed in the LTB 
model to confirm the quantitative agreement, but we do not expect a small spatial 
variation of the Hubble parameter to change the result notably. 

Indeed, the smallness (~ 15%) of the spatial variation in the Hubble parame- 



ter is another significant result of the fits in Sects. 34. and 3^. The variation can 
be considered small, as it is of the same order with the uncertainty of the model- 
independently 5 deduced value for the local Hubble rate [36]. The variation of the 
Hubble parameter found by Alnes, Amarzguioui and Gr0n [26] has also similar mag- 
nitude, but their model contains a large (~ 400%) variation in the matter density. 

In addition to the ones discussed in Sects. [3~1~1 - |3.5| , we have considered various 
other forms for the boundary condition functions VLm{t) and Hq{t). A generic out- 
come is that inhomogeneities in H (r) appear to have a much bigger effect on the 
goodness of the fit than the inhomogeneities in Qj^(r). In fact, we have been able to 
obtain good fits for the Hq(t) = H + AHe~ r ^ r ° -models with inhomogeneities of al- 
most any kind on ^(f). Moreover, we have found that models with Ho(r) = const, 
do not give a good fit irrespective of what kind of inhomogeneities we have inserted 
in f2jvf( r )- As an exception, a relatively high (O ~ 30) and thin (r ~ 150 Mpc) 
peak of the form 6 ^AffV) = fio e ~ r ^ r ° + ^o/(3e 4 ) seems to slightly improve the fit, 
giving x 2 ~ 180 when H (r) = const.; the result can be understood by looking at 
Eq. (|3.8| ), which places the derivatives of the functions VIm{t) and Hq{t) on a similar 
position, and as H' (r) < can mimic acceleration, so can Q' M (r) < 0. 

The results of Sect. |3]^ indicate that the inhomogeneities of the expansion rate 
and vacuum energy are mutually exclusive in the sense that their combination does 
not lead to a better fit. A good fit can be achieved either by having vacuum energy 
but no inhomogeneities; by having inhomogeneities but no vacuum energy; or by 
having both about the half amount of their separately deduced best fit values, i.e. 
they seem to have a very similar effect on the supernova observations. 

The early supernova data had a discrepancy between the Hubble parameter de- 



5 Note that the smaller uncertainties found in the CMB data analysis cannot be used here as 
those fits assume that the entire universe is perturbatively close to the homogeneous FRW model. 



3 The term f2o/(3e ) has been included to make the matter density in Eq. (B.£) positive 
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duced from the low redshift sample and the one deduced from the high redshift 
sample [37], but it seems to have vanished from the latest supernova data [38]. How- 
ever, the results of Sects. pTT| - pT5| suggest that the feature still exists, and rather than 
a discrepancy in the data, it is actually an alternative explanation to dark energy. 
Note that the LTB model certainly describes a universe with spatially varying ex- 
pansion rate more realistically than two separate FRW models, one for high redshift 
regime and another for low redshifts. In addition, the discrepancy in the early super- 
nova data between the local and global Hubble parameters had a similar magnitude 
with our best fit value for the variation of the Hubble parameter. The length scale 
associated with the variation of the Hubble parameter, found to lie within the range 
of about 200 to 1200 Mpc in our analysis, matches the explanation as well. 

As can be seen in Sect. |3.2| , the explicitly solvable LTB model with fij^r) = 1 
fits the data essentially as well as the more complicated ones with ^M( r ) 7^ 1- Hence 
it will be a useful example when discussing the time evolution of the inhomogeneities 



in Sect. [J77| and the average expansion in Sect. |. 

In order to respect the cosmological principle, we should not live in the dead 
center of a region that is expanding faster than the global average. It was shown in 
[39] that the supernova data does not impose severe restrictions for the location of 
an off-center observer, but as argued in [34], perhaps the most relevant constraint 
comes from the dipole anisotropy of the CMB. Here we give a rough estimate for the 
off-center distance allowed by the CMB dipole. 

We define the effective peculiar velocity, caused by the inhomogeneous expansion, 
that an observer at the coordinate r = d has relative to the symmetry center 

v p (d)=d(H (0)-H (d)), (3.6) 

which simply measures the deviation from the Hubble law. In a homogeneous uni- 
verse, the effective peculiar velocity of Eq. (|3.6| ) is identically zero and the observed 
peculiar velocity is accounted for motion relative to the comoving coordinates. We in- 
stead assume this coordinate velocity negligible and require that the effective peculiar 
velocity of Eq. ( |3.6| ) gives rise to the CMB dipole. Inserting the velocity deduced 
from the CMB dipole, v p = 370 km/s [40], and the boundary condition function 
H (r) of Eq. Q31|) with the best fit values r = 450 Mpc, AH = 10 km/s/Mpc to 
Eq. ( |3.6|) gives d = 140 Mpc, which is somewhat more than was found in a more 
thorough analysis for a different model 7 [34]. This is a non- negligible fraction of 
the size of the inhomogeneity, d/ro = 0.3, so according to this naive estimate the 
CMB dipole does not constrain our position strictly to the center. Of course even 
greater off-center distances are allowed if we would also have coordinate velocity in 
the opposite direction to counterbalance the effective peculiar velocity of Eq. (|3.6|). 



7 In [34] the function H (r) was fixed up to the distance of the LSS, whereas we have left its tail 



to be fixed from the CMB observations. Thus, instead of AH, we use H o (0) — H (d) in Eq. (3.6) 
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3.7 Time evolution of the inhomogeneities 

In the preceding discussion we defined the models by giving the boundary condi- 
tions Ho(r) and QmO") on the present-day spatial hypersurface. However, as is 
inherent in the Big Bang cosmology, the question of the "naturalness" of the ini- 
tial conditions arises in models with inhomogeneous expansion but a homogeneous 
present-day matter distribution. In this section, we address the problem by consid- 
ering the time evolution of the inhomogeneities in the explicitly solvable LTB model 
with QM{ r ) — 1- We expect the same qualitative features to hold in the more general 
case with flM{r) < 1 as well. For more on the subject, see e.g. [9]. 

Combining Eqs. ( |2.6| ) and ( |2.10| ) we obtain the following expression for the 
matter density: 



where 



Pm(t, t 



PAf(r,t) = p M (r,t ) 



3H2(r) n f v 



Al{r)A' {r) 
A 2 ( rj t)A'(r,t) 



With tt M (r) 



SnG 
1 and A (r) - 



A n (r) 



3A' (r 



Of, 



reduces to 



+ 2 



HJr) 



r, Eq. (3 
3# 2 (r) + 2rH' {r)H (r 



8nG[l 



2 



(f-f )][r(f-f„)fl&(r) 



3gp(r) 
2 



(*-*o))] 



(3.7) 



(3.J 



(3.9) 



which gives the explicit time dependence of the matter distribution. Using Eq. 
one finds that the ratio of the matter density near us (r = 0) to its asymptotic value 
(r = R > r ) is 

2 

(3.10) 



p M {0,t) 



pAi(R,t) 

At late times, t = T>t , Eq. ( 3.10|) gives 

Pm(0,T) 



H (0) (l + lH (R)(t-t )) 
H (R) (l + |jy (0)(t-*Q)) 



(3.11) 



p M (R,T) 

which means that, irrespective of the initial conditions, the homogeneity of the matter 
distribution is an attractor solution. 

By performing the integral of Eq. ( [2.14 ), one can calculate the explicit time 
dependence of the expansion rate: 

H (r) 



H(r,t) 



1 + 



, (t-t ) 
= 0) and global 



Eq. (|3.12|) gives the ratio of the local (r 

parameters 

H(0,t) H (0) (l + 'lH (R)(t-t )) 
H(R,t) H (R) (l + |^ (0)(t-t )) 



(3.12) 



R > r ) Hubble 



(3.13) 
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which tells us that the inhomogeneity of the expansion rate will also vanish in the 
late {t = T t ) universe: 

By comparing Eqs. ( |3.10| ) and ( |3.13| ) one sees that the inhomogeneity of the 
matter density will diminish more rapidly than the inhomogeneity of the expansion 
rate. This explains in a natural way why, with generic initial conditions, the present 
universe would have negligible inhomogeneities in the matter distribution but de- 
tectable inhomogeneities in the expansion rate, like our best fit model in Sect. VTA 



On the other hand, as noted in Sect. |3.6| , even if today's variations in the matter 
density were of the same order with the required small variation in the expansion 
rate, the observations might not be able to detect them. 

Turning the argument around, Eqs. (|3.10|) and ( |3.13 ) tell us that in order to 



have noticeable inhomogeneities today, the early universe should have been very 
inhomogeneous. Naturally, we do not have an explanation for this as there is no 
theory that would determine the initial conditions of the universe. However, when 
considering the backward time evolution, one should keep in mind that the dust 
approximation will break down at some stage and as a consequence, Eqs. (|3.9|) and 



( |3.12|) do not represent the inhomogeneities of the very early universe. 



3.8 The age of the universe 

The age of the perfectly homogeneous universe is a well-defined number, the differ- 
ence t age = t — tBB of the time coordinate today t and the time coordinate for which 
the scale function goes to zero a(t — > £_bb) 0. However, with inhomogeneities the 
notion of the age of the universe becomes more complicated as the age can depend on 
the spatial location (see e.g. [9]). In the LTB model this can be seen by noting that 
the singularity condition A(r, tsB) = now defines a curve tss(r) on the (r, £)-plane 
whereas in the homogeneous FRW case a(tBB) = defines a unique point tsB on the 
t-axis. We give the explicit expressions of the age of the LTB universe as a function 
of the boundary conditions H (r) and £7a/(V) for different cases in Appendix 

The expressions ( |A. 1 ) , ( JA.2|) , ( |A.3| ) and ( |A.4j ) are similar to the homogeneous 



FRW case, except for the r-dependence. In this case, the local value of the Hubble 
parameter gives the order of magnitude for the local age of the universe. Thus, if 
^m(^) = const., already small spatial fluctuations in the Hubble parameter H (r) 
lead to quite big variations in the age of the universe between different spatial lo- 
cations. For a constant ^AffV) the scale of the spatial variation in the age of the 
universe is simply given by (AH)- 1 , and at scales of 0(1000) Mpc would thus be 
some billions of years for the best fit models of Sects. |3.1| - |3.4j . Although theoretically 
inhomogeneous Big Bang could be a viable possibility (one need only to think of e.g. 
colliding branes; see [41,42]), age differences of such great magnitude could be hard 
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to reconcile with any existing model of the very early universe. They certainly con- 
tradict the spirit of inflation, where the whole observable universe arises from a single 
causal domain. However, the possibility remains that, for some reason, H Q (r) and 
Hjw(r) are related in such a fashion that they give t age (r) = constant. The inhomo- 
geneous LTB model would continue to fit the supernova data well as demonstrated 
by the model of Sect. 



4. Is accelerated average expansion needed? 

The expansion of the dust dominated LTB universe does not accelerate anywhere but 
nevertheless can, as shown e.g. in Sect. |3|, fit the supernova observations. However, it 
has been argued that a certain kind of average expansion can accelerate even though 
the local expansion would decelerate everywhere [13, 16, 18]. The aim in this section 
is to study if this kind of average acceleration appears in our best fit models of 
Sect. |3|. For this purpose, we need the Buchert equations that describe the averaged 
dynamics of a general irrotational dust universe [16] 

3^|| = -4nG( PM ) D + 2 -{{6 2 ) D - (9)1) - (a^) D (4.1) 

3^|| = 8nG(p M ) D - i<J2< s % - 1 -{(6*) D - (6)1) + \(a^a, v ) D (4.2) 
^(pm)d + S^(pm)d = 0, (4.3) 



where a^a^ > represents the shear, a{t) = (J A/det^j^x) 1 / 3 is the averaged 
scale factor, D is the integration domain, g^ is the spatial part of the metric, is 
the curvature scalar of the t = const, spatial hypersurface, 6 = V fl u tl is the expansion 
scalar and the spatial average of a scalar S is defined as 



J D S y/det\gij]d 



x 



{s)d s 7 /jrrr/3 • ^ 

The difference of Buchert's acceleration equation ( |4.1| ) and its homogeneous FRW 
counterpart is known as the backreaction 

Q = 1((0 2 )d- (^)-Kv> D . (4-5) 



The average expansion accelerates if the right hand side of Eq. (|4.1|) is pos- 
itive; this can be achieved by having large enough variance of the expansion rate 
although it is counterbalanced by the average shear. The variance gets large val- 
ues when contracting [9 < 0) and expanding (8 > 0) regions coexist, and in fact 
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the average acceleration has been connected to gravitational collapse [13,33]. How- 
ever, as demonstrated in [43], a globally expanding dust universe can have average 
acceleration as well. 

The definition of the shear tensor is 

0> = ~(V M u„ + V u u^) - ^(g^ + u^u u )V a u a , (4.6) 

where is the four- velocity of the dust. Calculating the shear and the expansion 
scalar for the LTB metric in the coordinates of Eq. ( |2.1| ) gives 

<""<*•=! (4 - 2 < 4 - 7 > 

9 = 2j + ^ = 2H(r,t) + H r (r,t) . (4.8) 

Using the field equation (|2.5|) one finds that the two Hubble functions H(r, t) and 
H r (r,t), defined in Eq. ( |4.8j ), are related at the present time 8 t as 

H r (r,t ) = H (r)+rH' (r) , (4.9) 

so that at t — to the shear and the expansion rate have the expressions 

o^o, u (t = t ) = ^{rH',{r)) 2 

9(t = t ) = 3H (r)+rH' (r) . (4.10) 



Inserting Eqs. ( |3.4| ) and ( |4.10| ) to Buchert's acceleration equation ( |4.1|) , we obtain 



the exact expression for the average acceleration of the model in Sect. 

= - 3 ^oH 2 + 6(H 2 ) + 4(rH' H ) - Q(H } 2 - A(rH' ) (H ) - 2 -(rH' ) 2 . (4.11) 
We choose the integration domain in Eq. ( 4.11 ) as the origin-centered ball with ra- 



dius R, taken to be greater than the inhomogeneity scale (e.g. R ~ 2r ) but smaller 
than the horizon distance R < 1/H. Note that at t = to we have a/ detf^] = 
r 2 / yT— k(r) in the integration measure. As the curvature function of Eq. (|2.11|) is 
small for the models of Sect. ^ we can make the linear approximation, 1 / y/l — k(r) ~ 
1 + k(r)/2, for the integrals in Eq. ( f4.4|) to get an approximate expression for the 
average of a scalar S 



(S)m<r) = Zr'T ?"'' - ~ ( S )c + 5(5*)c - ^(S)c(k)c , (4.12) 



J u Sr 2 (l- kir))- 1 / 2 ^ ^ wov l /ofN 1 
/ fl r 2 (l- kir))- 1 / 2 ^ 



8 In general, i = < is the moment when the boundary conditions Om(?") an d H Q (r) have been 
specified. 
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where we have used (S) c to denote the "coordinate" average: 

, . L Sr 2 dr 

(S) c = ^r—- • (4-13) 
J r z ar 

Using the approximation of Eq. ( 4.12j ), we can calculate the integrals of Eq. ( 4.11 ) 
to obtain an expression for the average acceleration of the model in Sect. 
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(4.14) 



2 24 288 x 1728 x 2 20736 x 3 

9 h 4 16 h 3 9 
"256^3 + 81^3 + 2^ ^'x 3 
where we have denoted h = AH/H and x = R/r . A careful inspection reveals that 
the negative first term on the right hand side of Eq. ( [4. 14 ) dominates, at least for 
the values of the parameters r , AH, H and flo that give a good fit. Thus, according 
to Eq. ( 4.14j) , the average expansion does not accelerate in the best fit model of 
Sect. |3.4| . The validity of Eq. ( 4.14| ) is verified by the fact that it actually deviates 
only about 10~ 7 from the numerically computed accurate result of the integrals in 
Eq. (|4~TT| ) for the best fit values h = 10/57, r = 450 Mpc and f2 = 0.4. We have 



repeated the computation of the average expansion also for the model of Sect. [3J. 
and found no acceleration there either. 

Finally, let us consider the LTB model with VLm{t) = 1. For this model, both 
the backreaction of Eq. ( |4.5| ) and the 3-curvature vanish identically [43]. Thus, 
the Buchert equations ( [4.1|) , ( |4.2| ) and (|4.3|) reduce to the FRW equations of the 
flat and homogeneous dust universe. However, a creature living in the LTB universe 
would observe the luminosity-redshift relation calculated from the exact Einstein 
equations, not from the FRW equations. The difference between these results is 



significant since, as shown in Sect. |37J, the LTB model with fijw(r) = 1 fits the 



supernova observations whereas the flat matter dominated FRW model does not [2]. 
As the averaged matter density (pm) d has only a small dependence on the averaging 
domain D in our models, the conclusion would not change even if one were to use 
a different smoothing scale Ri for every supernova at r — R4. Overall, there seems 
to be no meaningful way of using the averaged equations to calculate the actual 
observables, e.g. the luminosity-redshift relation. 
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5. Conclusions 

We have studied the effect of inhomogeneities on the cosmological supernova ob- 
servations in the spherically symmetric matter dominated LTB universe. Our main 
conclusions follows: 

1. The inhomogeneous matter dominated LTB model fits the supernova data bet- 
ter than the homogeneous FRW model with the cosmological constant. 

2. Inhomogeneities of the expansion rate seem to have much bigger effect on the 
goodness of the fit than the inhomogeneities of the matter distribution. 

3. The model giving the best fit has perfectly uniform present-day matter density 
with Qm ~ 0.3, in agreement with the galaxy surveys, but its Hubble parameter 
varies from the local value of ~ 70 km/s/Mpc to the value of ~ 60 km/s/Mpc 
over the distance scale of ~ 500 Mpc from us. 

4. Neither local nor average acceleration is needed for a good fit to the supernova 
data. 

5. At least in the LTB universe, the averaged Buchert equations give a notably dif- 
ferent prediction for the observable luminosity-redshift relation than the exact 
Einstein equations. 

The first of the above listed results is neither new nor surprising, as it was proved 
in [14] that any isotropic set of observations can be explained by appropriate inho- 
mogeneities, and later demonstrated with examples by several authors. Instead, the 
second result is a new one, as is the model that explains both the observations of 
supernovae and the observed homogeneity in the large scale distribution of galax- 
ies without invoking dark energy. Moreover, we found in Sect. [3.7| that within the 
dynamics of the Einstein equation, the universe may evolve to this kind of a config- 
uration without any substantial amount of fine-tuning in the initial conditions. In 
addition, as the supernovae and galaxy surveys constrain the boundary condition 
functions H (r) and VLm{t) only for the small values of the radial coordinate r, the 
tails of these functions can still be freely chosen to give a good fit for the CMB power 
spectrum as well. In fact, Alnes, Amarzguioui and Gr0n have realized this in practice 
by fitting the LTB model to the first peak of the CMB [26]. It has also been shown 
that even a homogeneous model can fit the CMB data without dark energy [44]. 

Our another focus was on the averaged dynamics of the dust universe, given 
by the Buchert equations that illustratively demonstrate the potential cosmological 
significance of the inhomogeneities. In Sect. [|, we showed that none of our dust 
models that fit the supernova observations have accelerated average expansion. On 
the other hand, one can easily construct examples that have average acceleration 



- 22 - 



but are ruled out by observations, whereas dust models that would both have the 
average acceleration and fit the cosmological data are still to be presented. Moreover, 
as discussed at the end of Sect. f|, it seems that the Buchert equations cannot be 
used to calculate the observable luminosity-redshift relation. Although explicitly 
demonstrated only in the LTB universe, we expect that it might be the case in 
the more general models as well; after all, the Buchert equations are obtained by 
integrating over a spatial hypersurface whereas the observables are calculated by 
integrating over the past light cone. Altogether, averaging over a spatial hypersurface 
is perhaps not practical if one wants to make comparison with the observations. 
Instead, probably a better alternative is to relate the inhomogeneities directly to the 
observations, as done e.g. in this paper and, if needed, take the averages only of the 
relevant observables, such as the energy flux and redshift of light. 

Of course, some sort of smoothing is implicitly assumed also in the LTB model. 
Therefore, the study of observations in the LTB universe does not answer the in- 
teresting question about the potential cosmological consequences of the small scale 
lumpiness caused by galaxies. Their effect could be significant, possibly depending 
also on how smoothly dark matter, perhaps the dominant energy component of the 
universe, is distributed in the universe. 
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A. The age of the LTB universe 



The age of the LTB-universe can be calculated by integrating the field equation 
( |2.12|) . In the cases where the result is an elementary function of the boundary 
conditions, the results are: 

1. tt M (r) < 1 and fi A (r) = 0: 



, l-lMQ-^arsinh^^ 

2. Q M (r) > 1 and fi A (r) = 0: 

fi A/ (r)arcsin v / ^^i - y/Q^Fj - 1 
H (r)(tt M (r)-l)^ 



t age (r) = T^T, • (A.2) 



3. fl M (r) = 1 and Q A (r) = 0: 



= — . (A.3) 



4. n M (r) + n A (r) = 1: 



4ge r = oh- / \ /, = r=s ■ A - 4 
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